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Abstract. 

We report on our recent evaluation of the two-photon width of the pseudoscalar quarko- 
nia, r] c (nS) and rjb(nS) in an approach based on Heavy-Quark Spin Symmetry (HQSS). To 
what concerns the IS state rj Cl our parameter-free computation agrees with experiments, 
as well as most of other theoretical works. On the other hand, our computation for the 
25-state looks 2S like a confirmation that there may exist an anomaly related to the decay 
of rj' c , especially in the light of the new preliminary result of the Belle collaboration. We 
also point out that the essentially model-independent ratio of r/b two-photon width to the 
T leptonic width and the r\\, two-photon width could be used to extract the strong coupling 
constant a s . 
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1. INTRODUCTION 

Since the J/ ip discovery, bound states of heavy quarks -the heavy quarkonia- are 
expected to provide physicists with ideal means to study the main properties of 
QCD. With the time though, it appeared progressively that such quark systems are 
not so easy to understand and controversies about their production mechanisms 
are still going on (TJ [2] . Fortunately, the physics involved in their decay seems to 
be rather well understood within the conventional framework of QCD [TJ [3]. 

However, recently two experimental estimations of the two-photon width of 
the r]' c , one published by the CLEO collaboration jl] (T 7J (r]' c ) = 1.3 ±0.6 keV) 
and another preliminary by the Belle collaboration (r 77 (7^,) = 0.59 ± 0.13 ± 0.14 
keV) [5j |6] contradict most of the existing theoretical predictions lying in the range 
3.7 to 5.7 keV El El E2 ED E3 E2 E]. This is rather surprising since such 
electromagnetic decays of non-relativistic states should be rather easy to describe 
from a theoretical point of view. 

It was therefore our purpose in [T5J [16] to reanalyse such a process from a very 
basic starting point: the heavy-quark spin symmetry (HQSS). Indeed, in a non- 
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relativistic system, the difference between the IS and 2S widths would appear only 
at the level of the wave function at the origin. In virtue of such a symmetry, both 
wave functions should not differ much from the ones of 3 Si state, which are a priori 
well known from the leptonic decays of J/ip and ip' . 

We report here on our effective approach based on HQSS and compare its 
results with other available calculations. Since none of them is able to predict 
both experimental measurements for r 77 (r] c ) and T^(r]' c ), we also discuss some of 
the hypotheses made by CLEO and Belle for the extraction of those widths from 
their experimental observables. We finally report on the corresponding predictions 
for rjb(nS) states. 

2. EFFECTIVE LAGRANGIAN FOR % DECAY INTO 

TWO PHOTONS 

In the two-photon decay of a heavy quarkonium bound state, the outgoing-photon 
momentum is large compared to the relative momentum of the quark- ant iquark 
bound state, which is 0(A/m,Q), with A <C uiq. One obtains an effective Lagrangian 
for the process QQ — ► 77 (represented by the first diagram in FigJlJ) by expanding 
the heavy-quark propagator in powers of q 2 /mQ, and neglecting 0(q 2 /rriq) terms 
(q = Pq ~Pq)- Like leptonic decay, the two-photon decay, in this approximation, is 
described by the following effective Lagrangian: 

Cf a = -c 2 (QYQ)(V) (1) 

with 




FIGURE 1. Effective coupling between a QQ and two photons (left) and a lepton pair (right) 

The factor 1/(M? Sq + bis Q Mi So ) in c\ contains the binding-energy effects (the 
binding-energy b is defined as b = 2mq — M) and is obtained from the denominator 



of the heavy-quark propagator (k\,k2 being the outgoing-photon momenta): 

1 



[{k 1 -k 2 ) 2 /A-ml 



(3) 



The decay amplitude is then given by the matrix element of the axial- vector current 
Ql^l^Q similar to the quarkonium leptonic decay amplitude which is given by the 
vector-current matrix element Q^Q for 3 Si — > £ + £~ . Thus for decays of S"-wave 
quarkonium into two photons or a dilepton pair ££, we have: 



Ms Sl 

M„ = -MQ 2 Q (Ana em ) £* e^egeffiff (4) 

M 1So + bi So lVh So 



where 



(0\QYQ\ 3 Si) = f 3Sl Ms Sl e^ (0\QYl 5 Q\ 1 S )=tfi So P". (5) 
from which the decay rates are: 



1 U\ *->lJ — TTVlr ) 1 77V ^O) ~ jT, • 



3M 35l -ta-u, MiSo 



(6) 



By taking Mss 1 f 2 Sl = 12|^(0)| 2 , we recover the usual non-relativistic expression for 
the decay rate which, whose NLO QCD radiative corrections are given by 

r NLO ( 3 s 1 ) = r LO ( 3 s 1 ) (1-^) (7) 

V 7T 3 / 



a s (20 



7T 



2^ 



3. MATRIX ELEMENTS OF LOCAL OPERATORS 

We have shown that in the approximation of neglecting 0(q 2 /rriq) terms, the two- 
photon decay amplitude is given by the 1 S'o decay constant fis - We now derive a 
symmetry relation between fi So and f3 Sl , the 3 Si leptonic decay constant using the 
relativistic spin projection operators for a relativistic Bethe-Salpeter quarkonium 
bound state [T7] . 

Consider now the matrix elements of local operators in a fermion-antifermion 
system with a given spin S and orbital angular momentum L 



A = J ^TrO(0) X (P,q) (9) 

P is the total 4-momentum of the quarkonium system, q is the relative 4-momentum 
between the quark and anti-quark and x{PiQ) ls the Bethe-Salpeter wave function. 



For a quarkonium system in a fixed total, orbital and spin angular momentum, 
x{P,q) is given by (q is the relative 3-momentum vector of q) 

X (P,q;J,J z ,L,S) = 2n5(q°-^)^ LM (q)(LM;SS z \JJ z ) 



:J 3 E <P/2 + q, s)v{P/2 - q, s) <~s; ~s\SS,) 
V 2 2 

= 27r8(q -^-)^ LM (q)V S s z (P,q){LM;SS z \JJ z ) . (10) 

M,S Z 1 m 

The spin projection operators Vss z (P,q) are 



^o,o (P, q) = \ ~ i~( p / 2 + m hs K/72 " rf) + m] 



5m 3 



*W^) = \ riH/'/H^ + m^K^-^ + m] • (11) 

V 8771 

For 5-wave quarkonium in a spin singlet 5 = and spin triplet 5 = 1 state: 

(2tt) 

In this expression the g-dependence in the spin projection operator has been 
dropped and the integral in Eq.(12) is the 5- wave function at the origin [T9] : 

<] '' ^oo(g) = -^=n (o) . (13) 



A( 2S+1 Sj) = Tt(O(0)V JJz (P,0)) J T^lMff) • (12) 



Using Eq.(ll) and Eq.(12) to compute the matrix elements {^\Q^^Q\P) and 
(Q\QlnQ\V) for the singlet 5 = pseudo-scalar meson P and for the triplet 5 = 1 
vector meson V, we find, neglecting quadratic 0(q 2 ) terms. 



Ws^' '' 4 ™'- Wh^^t^ (14) 

Thus we get the relation 

f 1S0 ^f 3Sl +O(b 2 /M 2 ). (15) 

It is expected that this relation holds also for excited state of charmonium and 
bottomonium where the binding terms 0(b 2 /M 2 ) can be neglected. This is a 
manifestation of heavy-quark spin symmetry(HQSS). In this limit, the two-photon 
width of singlet 5 = quarkonium state can be obtained from the leptonic width 
of triplet 5 = quarkonium state without using a bound state description. This 
approach differs from the traditional non-relativistic bound state approach in the 



use of local operators for which the matrix elements could be measured or extracted 
from physical quantities, or computed from QCD sum rules [20|. l2Tj and lattice QCD 
[22]. 

The ratio of the rj c two-photon width to J/ip leptonic width in the limit of HQSS 
is then: 

^^ 3Ql ^( 1 + ^I^>). (16) 

IW/VO M ric \ 7T 3 J 

For M Vc = Mj/ij,, the above expression becomes the usual non-relativistic result 
[231 E3] as mentioned above. From the measured J ftp leptonic width, we get 
F-yyiVc) = 7.46 keV. Including NLO QCD radiative corrections with a s = 0.26, we 
find r 77 (?7 c ) = 9.66 keV in near agreement with the world average value 7.4 ± 0.9 ± 
2.1 keV. A similar result is obtained in [23] which gives 8.16 ±0.57 ±0.04 keV . 

Thus the effective Lagrangian approach successfully predicts the rj c two-photon 
width in a simple, essentially model-independent manner. 



4. HQSS PREDICTIONS FOR T^rQ 

To obtain the prediction for rf c , we shall apply HQSS to 2S states. Thus, assuming 

f 2 , 

f'v'c = fi>'i an d neglecting binding-energy terms, we find: F (rf c ) = F (rj c )j^- = 

3.45 keV. This value is more than twice the evaluation by CLEO and five times 
the one by Belle. Our results is however nearly in agreement with other theoretical 
calculations [UEIE] as shown in Table 1. Other approaches [TQl EH E21 EH] seem to 
be closer to the latter measurements but then undershoot clearly the measurements 
for 77 c . 

Including binding-energy terms, for M Vc ~ Mj/^, M v ' c ~ ilfy, we have 
which gives 

r 77 (77 c ) = 4.1keV . (18) 



TABLE 1. Theoretical predictions for T 17 (t] c ) and r 77 (r^). (All values are in units of 
keV). 
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11.8±0.8±0.6 
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5.5 6.2 




3.5-4.5 


3.7 


4.44 ±0.48 


5.7±0.5±0.6 


1.38±0.3 


2.1 


1.8 3.36-1.95 



Binding-energy corrections seem then to worsen comparison with data and 
maybe point at anomaly in the decay r) c (nS) —>■ 77. However, before drawing such 
a conclusion it is necessary to discuss the experimental hypotheses made for the 
extraction of the aforementioned widths. 



5. B(r] c (nS)^KKir) 
The measured value from CLEO [I] : 

r 77 (^) = 1.3±0.6keV, (19) 
was effectively done by considering the following quantity : 

RK/Vc) = ' 77 ^! X !!^^ T ! = °-18 ± 0-05 ± 0.02. (20) 

To obtain r 77 (?■/£) from the above data, they made the assumption that 

B( V ' C -> K s Ktt) « B( Vc -> ^Kvr) (21) 



and in turn found the result of Eq. (19). 

Such an assumption is in fact supported by a couple of observations. Belle 
measurements of B — > r] c K and B — > gives 



= B(B^Ik^B(!I^KsK + «-) = 0-38 ±0-12 ±0.05. (22) 



Using the approximate equality Eq. (21), one would obtain 

gwg sG4 (23 ) 

which agrees more or less with the QCD factorization (QCDF) prediction J27J : 

On the other hand, it is expected from Sf/ (2) flavor symmetry that one would 
have the approximate equality between the ratios 

B(B° -> rfJC ) „ B{B+ -> 

B(B° ^ Vc K°)~ B(B+ ^r] c K+)' 1 J 

This is indeed the case since the following ratio obtained from BABAR 

B{B+ -> 77^+ 



0.38 ±0.25 (26) 



corresponds again to Eq. ( 23 ) . 

Those observations therefore tend to support the assumption of the approximate 
equality between the rf c — > KKtt and r] c —>■ KKn branching ratio. This would 
confirm a small r]' c —>■ 77 decay rate as quoted above. 



Albeit, we also note that the good agreement with QCDF predictions for the 
measured ratio B{B° -> r]' c K°)/B{B -> Vc K°) and B{B + -> V ' c K + )/B{B + -> r/ c K+) 
at Belle and BABAR suggests that f v > c /f Vc ~ /<// / /,//</>, which in turn supports HQSS 
and our predicted value for the two-photon width which is more than twice bigger 
than the CLEO estimated value shown above. More precisely, comparing R(r]' c /r] c ) 



with R(j]' c K J T] C K) and using QCDF value given in Eq.(24), we find 

R(rfJr h )»R(rf e K/r h FC)/0.9 . (27) 
The Belle data in Eq. d22b would then implies R(v' c /Vc) « 0.42 ±0.13 ±0.05, twice 



bigger than the CLEO data shown in Eq. (20). 

Since QCD gives B(r) c — > 77) ~ ^(v'c ~^ 77jand the predicted B{r\ c — > 77) agrees 
well with experiments, we expected that the large value for the measured r)' c total 
width would imply a large value for the rf c two-photon width. Thus it is difficult 
to understand the very small recent Belle measured rj c two-photon width. 



6. HQSS PREDICTIONS FOR r 77 (r/ 6 ) AND V^(rf b ) 

Since the 6-quark mass is significantly higher than the c-quark mass, the effective 
Lagrangian and HQSS approach should work better for bottomonia decays to 
leptons and photons. We thus have: 

(neglecting the small b Vb /M Vb binding-energy term). This gives T^(r] b ) = 560 eV 
(a s (M T ) = 0.16, M Vb = 9300 MeV) . 

For r]' b and higher excited state, one has (M Vb ~ M T and M v > ~ M T '): 

which gives r 77 (?7^) = 250eV and r 77 (r/^) = 187eV. In Table. [2] we give our pre- 
diction for the two-photon width of r] b , rf b and r]' b ' together with other theoretical 
predictions. We note that our predicted values are somewhat higher than other 
predicted values. 

Eq. ( 28 ) can be used to determine in a reliable way the value of a s . The momen- 
tum scale at which a s is to be evaluated here could be in principle be fixed with 

Rr) b - 

Further check of consistency of the value for a s may be possible in future 
measurements on the rjb and its two-photon decay branching ratio: 

^ = ?<S%l-7.8*). (30) 
T gg (r] b ) 2 aj \ n J 



TABLE 2. Summary of theoretical predictions for r 77 (r7b), T 17 (ri' b ) and r T0 , (77") . (All 
values are in units of cV) . 
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7. CONCLUSION 

We have shown here that effective Lagrangian approach and HQSS can be used to 
compute quarkonium decays into leptons and photons with relativistic kinematics, 
for both ground states and excited states of heavy-quarkonium systems. 

We emphasised with our basically model-independent calculations that either 
HQSS holds for radially excited charmonia and there is a not-yet-understood 
specificity in the decay r( c — * 77; either HQSS is strongly broken and this hints 
at large relativistic corrections at work in such decays . This would in turn explain 
why most of available models on the market are unable to give correct predictions 
for both T] c and r]' c decays. 

Measurements of the two-photon widths for rjb and higher excited states could 
provide with a test for HQSS and a determination of the strong coupling constant 
a s at a scale around the T mass, similarly to what has been done with the T 
leptonic width in the past. 
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